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We defend a natural division of the energy density, energy flux and momentum density of elec- 
tromagnetic waves in linear media in electromagnetic and material parts. In this division, the 
electromagnetic part of these quantities have the same form as in vacuum when written in terms 
of the macroscopic electric and magnetic fields, the material momentum is calculated directly from 
the Lorentz force that acts on the charges of the medium, the material energy is the sum of the 
kinetic and potential energies of the charges of the medium and the material energy flux results from 
the interaction of the electric field with the magnetized medium. We present reasonable models for 
linear dispersive non-absorptive dielectric and magnetic media that agree with this division. We 
also argue that the electromagnetic momentum of our division can be associated with the electro- 
magnetic relativistic momentum, inspired on the recent work of Barnett [Phys. Rev. Lett. 104, 
070401 (2010)] that showed that the Abraham momentum is associated with the kinetic momentum 
and the Minkowski momentum is associated with the canonical momentum. 
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I. INTRODUCTION 

The discussion of the correct form of the energy- 
momentum tensor of electromagnetic waves in material 
media has been the subject of a hot debate for more 
than one century [1]. The two most famous proposals 
were made by Minkowski and Abraham on the beginning 
of last century ^ and have distinct predictions for the 
momentum of light in a medium. The Abraham formu- 
lation predicts that an electromagnetic wave with total 
momentum p in vacuum decreases its momentum to p/n 
after entering a dielectric medium with refraction index 
n, while the Minkowski formulation predicts that the mo- 
mentum of the wave is increased to np after entering the 
medium. Several experiments on this subject were made 
in the past 60 years [5l411| , many of them with the pur- 
pose of solving the debate, but the fact is that no decision 
can be made on experimental grounds. As it was pointed 
out by Penfield and Haus [T^, the Abraham, Minkowski 
and other electromagnetic energy-momentum tensors are 
accompanied by distinct material energy-momentum ten- 
sors. And when the properties of the material energy- 
momentum tensors are carefully taken into account, both 
formulations predict the same experimental results. A re- 
cent review of this debate can be seen in Ref. [1]. The 
eventual conclusion is that there are many compatible 
ways of dividing the total energy-momentum tensor of 
the system into an electromagnetic tensor and a material 
one. 

Nevertheless, we can ask whether there is, among all 
possible ways of dividing the momentum and the energy 
of electromagnetic waves in material media into electro- 
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magnetic and material parts, a division that could be 
considered as more natural. Recently, many attempts 
have been made to calculate the momentum transfer from 
electromagnetic waves to material media directly from 
the Lorentz force law [I3H24] . We believe that this is the 
most natural way of calculating the material momentum. 
In Ref. [23] it is presented a review of some treatments 
based on Lorentz force calculations for experiments on 
the topic, showing that consistent results are obtained. 
We also believe that the most natural form for the mate- 
rial energy density is the sum of the potential and kinetic 
energy densities of the charges of the medium, which can 
be modified by the presence of an electromagnetic wave. 
The goal of this paper is to show, using reasonable mod- 
els for material media, that if the material parts of the 
momentum density and of the energy density of electro- 
magnetic waves in linear media are to be calculated in 
this way, the electromagnetic parts of the energy den- 
sity Uc.m., of the energy flux Sc.m. and of the momentum 
density Pc.m. must be 

Ue.m. = ?|EP -I- -^|B|2 , Se.n, = — E X B , (1) 

Pc.m. = eoE X B , (2) 

where E represents the macroscopic electric field, B the 
macroscopic magnetic field, Eq the permittivity of free 
space and fiQ the permeability of free space. These forms 
are identical to the expressions in vacuum, but written in 
terms of the macroscopic fields. The form for the electro- 
magnetic momentum density above had appeared before, 
for instance in the book of Livens [25,,. However, the dis- 
cussion of the energy of electromagnetic waves in linear 
media is practically absent from the literature. 

In Sec. [n] our proposal for the division is more care- 
fully described and compared with previous treatments. 
In Sees. |III[|IV| and[V| different models for dispersive non- 
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absorptive dielectric and magnetic media are shown to be 
consistent with the proposed material energy density. In 
Sec. |Vl]we reproduce the model introduced by Penfield 
and Haus j26| for the material energy flux and material 
momentum of an object with magnetic dipole moment 
in an applied electric field, showing the consistency of 
the division. This momentum of a magnetic dipole in 
an electric field is frequently called "hidden momentum" 
[371 and is a relativistic effect. In recent papers, Bar- 
nett et al. pOl [29l [30] associated the Minkowski mo- 
mentum with the electromagnetic canonical momentum 
and the Abraham momentum with the electromagnetic 
kinetic momentum. Based on the model of Penfield and 
Haus, we suggest that the momentum density Pc.m. of 
Eq. ([2]) should be associated with the electromagnetic 
relativistic momentum. Finally, in Sec. VII we present 
our concluding remarks. 



II. DISCUSSION OF THE DIVISION 

The microscopic Lorentz force density that acts on a 
density charge p' moving with velocity v, forming a cur- 
rent density j' = p'v, is 



f = ep' +j' X b 



(3) 



where e and b are the microscopic electric and mag- 
netic fields. If the separation between the charges in the 
medium is much smaller than the wavelength of the in- 
cident electromagnetic wave, we can associate averaged 
electric and magnetic fields E and B to volumes contain- 
ing a huge number of particles, but still with dimensions 
much smaller than the radiation wavelength. The mi- 
croscopic electric and magnetic fields can be written as 
e = E + ^e, b = B + where 5^ and (5b depend on 
the microscopic distribution of charges and average to 
zero. So the averaged force density in this volume can be 
written as 



f = Ep + J X B 



(4) 



where p represents the averaged charge density and J 
the averaged current density in the volume. In a lin- 
ear medium without free charges, we can write J — 
dP/dt + \/x M and p = -V-P, where P and M are the 
polarization and magnetization of the medium |31| . The 
electric and magnetic macroscopic fields are the fields 
that appear on the Lorentz force equation above, so these 
are the fields that can transfer momentum and energy to 
material macroscopic bodies. From our point of view, 
these should be the fields that transport electromagnetic 
energy and momentum. For this reason, these are the 
fields that appear in Eqs. ([T]) and ^ for the electromag- 
netic energy density, energy flux and momentum density. 
The Abraham and Minkowski momentum densities, on 
the other hand, are pAbr = E x H/c^ and PMin = D x B, 
respectively, where H = B//io — M, D = EqE + P and 
c = 1/ y/fM)£o is the speed of light in vacuum [2 [5] . 



Some of the previously cited papers that calculate the 
material momentum by the Lorentz force consider that 
the electromagnetic momentum density should have the 
Abraham value E x H/c^ [H [T71 dll Hg. Our pro- 
posal for the electromagnetic momentum density of Eq. 
([2]) agrees with the Abraham form in a non-magnetic 
medium. As it was discussed by Mansuripur |5T1[12], to 
have momentum conservation considering the Abraham 
form for the electromagnetic momentum density in mag- 
netic media, a modification must be done in the Lorentz 
force law. Mansuripur opts for this modification in the 
Lorentz force to avoid the concept of "hidden momen- 
tum" , that we will discuss in Sec. |VI| In a recent paper 
|23) . we showed that the form of Eq. ([2| for the elec- 
tromagnetic momentum density and the form of Eq. Q 
for the Lorentz force law are compatible with momentum 
conservation in a series of examples [36] . We also showed 
compatibility with the Balazs gedanken experiment |32j , 
showing compatibility with the theory of relativity when 
taking into account the "hidden momentum". Because 
in our formulation there is no need for a modification in 
the Lorentz force law, we believe that our treatment is 
more natural. 

Now let us discuss the energy of electromagnetic waves 
in linear media. In a closed system, the local conservation 
of energy can be written as 



duto 
~dt 



= -V • St 



(5) 



where wtot and Stot represent the total energy density 
and total energy flux of the system. For a system com- 
posed by an approximately monochromatic electromag- 
netic wave with frequency ujc in a linear dispersive non- 
absorptive dielectric and magnetic medium, we can write 



(Wtot) 



d{uJXc 



duj 

djoJXm) 
did 



2 

Mo(|Hp> 



E X H 



(6) 
(7) 



where (A) represents the average of A over one period of 
oscillation of the field, Xc{^) is the electric susceptibility 
and Xni(aj) is the magnetic susceptibility of the medium. 
So, if the electromagnetic parts of the energy density and 
of the energy fiux of electromagnetic waves in linear me- 
dia are given by Ec^. ([T]), the material counterparts must 
be given by 



(Umat) = 



Xo + W 



dXc 
duj 



da„ 



duj 



(|B| 



2mo 



-E X M 



(8) 
(9) 



with = Xm/(1 + Xm), SO that (Uo.m. + M„iat) = ("tot) 

and Sc.m. + Sniat = Stot- In the next three sections we 
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show that commonly used models for dielectric and mag- 
netic linear media acquire the energy density of Eq. ([s]) 
with the presence of an electromagnetic wave, and in Sec. 
VI we reproduce the model of Penfield and Haus [53] that 
justifies the material energy flux of Eq. ^. 



III. CLASSICAL MODEL FOR A LINEAR 
DIELECTRIC MEDIUM 

The first model that we present is the Drude-Lorentz 
model for dielectric media. A similar treatment was pre- 
viously made by Haus and Kogelnik [53]. In this model, 
the atomic nuclei are fixed and the electrons are bound 
by harmonic potentials Ui{ri) = mujfrf/2, where m rep- 
resents the electron mass, the distance of each electron 
to its equilibrium position and uji the natural frequency 
of the harmonic oscillator. Considering a non-absorptive 
medium, the equation of motion of the x component of 
each electron can be written as 



cPxi 



qE cos{ujt) 



m 

dxi 



(10) 



where q represents the electron charge and E — 
E cos{ujt)x is the incident electric field. The solutions 
of the above differential equation are 



q 



-E cos(aji) 



(11) 



Because we are disregarding the absorption of radiation, 
we must consider that uj is far from uji. 

If we call Ni the density of electrons subjected to the 
potential Ui, the electric susceptibility of the medium can 
be written as 



Xe(w) 



lEl 



eqE coa{ujt) 



E 



(12) 

The material energy density corresponds tho the sum 
of the kinetic (mi;f/2) and potential {mujfxf/2) energies 
of the electrons multiplied by their densities Ni. Using 
Eq. ( 11 ), we have 



^mat 



(13) 



Using Eq. (12), it is straightforward to see that Eq. 
(131 can be written as the first term of Eq. (|8|: 

Xc 



^mat 



- W 



(14) 



IV. 



QUANTUM MODEL FOR A LINEAR 
DIELECTRIC MEDIUM 



Now let us consider a quantum model for a linear 
non-absorptive dielectric medium. The medium is com- 
posed by molecules with fixed nuclei and Z electrons with 



charge q and mass m around. The Hamiltonian that cou- 
ples these nuclei and electrons may be very complex, but 
our results do not depend on its specific form. The eigen- 
states of the molecular Hamiltonian will be denoted !$„), 
and the respective eigenvalues /iw„. An arbitrary pure 
state of the molecular system can be written as 



|$(t))^^c„(t)|$„)e— 



(15) 



with EJcn(t)P-l. 

An electric field E — E cos{ujt)x acts as a perturbation 
on the system. Our treatment is based on Loudon's book 
[34) . The total Hamiltonian can be written as 

H = Ho + H'{t) , with i?o|*„) = /iw„|$„), (16) 

H'{t) = -Ecos{ojt)ZqXcm, (17) 

where icm is the x component of the center of mass po- 
sition operator of the Z electrons and Zq is the total 
charge of the electrons. 

From the Schrodinger equation with the Hamiltonian 
above applied on the state (151, after taking the scalar 
product on both sides with e*'^'"*($m|, we obtain 



Cmit) 



J2cn{t){<Pm\H'\'fn)e' 



(18) 



Let us consider that the molecule is in its fundamen- 
tal state |<I>o) before the interaction with the field starts 
and that the field produces a small perturbation, such 
that Co ~ 1 always. So we will consider only the term 
n = in the summation above to the calculus of Cm{t)- 
The electric susceptibility refers to the situation where 
the interaction between the field and the molecules has 
achieved an equilibrium condition, so we can take an in- 
definite integral as the solution of Eq. (18): 



Cra{t) 



EZqXjnO 



-j't(cj„i ^i{ujjn —uj)t 



UJ 



(19) 



with Xrao = ($m|^cm|$o), ^0 = 0, cq ~ 1. Wc uscd the 
definition Xqq = to define the origin of the coordinates 
system. 

The oscillating dipole moment of the molecule can be 
written as 



d^{t) = ($(t)|Z<7iem|$(i)) 

y^c^{t)ZqX*^^e~^--' + 



c.c. 



(20) 



where we disregarded terms CmC„ with m and n both 
different from zero, cq was considered 1 and c.c. stands 
for the complex conjugate. 

Using Eqs. (20 1 and (19), we can see that a medium 



composed by a density N of these molecules has an elec- 
tric susceptibility 



Xo(w) = 



£o|E| 



E 



X. 



mO\ 



1 



1 



(21) 
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Because we are considering a non-absorptive medium, 
Xe(w) is real. This approximation of non-absorptive me- 
dia is valid only when uj is far from all a;,,„ for which 
XmO 7^ 0. 

The material energy density is the expectation value 
of the non-perturbed Hamiltonian multiplied by the 
density of molecules, because the Hamiltonian con- 
tains the kinetic and potential energy of the electrons: 



^mat 



(22) 



Using Eq. 19 we have 



'm4 



mO| 



[w^ cos^(a;i) -|- sin^(a;t)] 



(23) 

Taking the temporal average over one period of oscil- 
lation of the field we have 



2^2 



(24) 



Using Eq. (21 1, it is straightforward to show that the 



energy density can be written as the first term of Eq. 



(Mmat) 



Xc +0; 



dco 



(25) 



If we have different kinds of molecules composing the 
medium, each one identified by a symbol (j), the electric 
susceptibility and the energy density of the medium will 
be 



Tot 



Xo 



,Tot\ 



(26) 



and Eq. ( 25 I remains valid. 



V. 



QUANTUM MODEL FOR A LINEAR 
MAGNETIC MEDIUM 



The quantum model that we use for a magnetic linear 
medium is based on Van Vleck's book [3S] and considers 
the response of electrons under the action of a central 
potential to an applied magnetic field. We use cylin- 
drical coordinates (z,s, 0) in this section. To simplify 
the treatment, we will consider that the magnetic field is 
generated by a vector potential A = Bs4>/2, such that 
B = V X A = _Bz, instead of considering the vector po- 
tential of an electromagnetic wave. The Lagrangean of an 
electron of mass m and charge q under the action of this 
magnetic field and an atomic central potential 3>(z, s) is 



C = ^m|rp + (/A • r — (/<I>(z, s) 



mz 



ms 



ms 4) 



2a2 



qs^ 



g$(z,s)(27) 



The Lagrangean does not depend on the variable 0, so 
the momentum conjugated to this coordinate is a con- 
stant of the movement of the system: 



— — r — ms (p 



qs'^B 



= constant. 



(28) 



So the z component of the angular momentum of the 
electron, = ms'^(j>, is a function of the applied magnetic 
field: 



qs'^B 



(29) 



where Iz (0) represents the original angular momentum of 
the electron, when the magnetic field is zero. 

The ratio between the dipole moment and the orbital 
angular momentum of an electron is q/{2m) |31j . so 
the variation of the electron dipole moment is A/i^ = 
^i.{B) - Ai.(O) = -q''{s^)B/im, with = 
I ^P) being the unperturbed quantum state of the electron. 
Considering that the medium has atoms with random 
orientations, on taking the average of the orientations 
we can consider, for each quantum sate of the electrons, 
Ejihmj = and = E,2/3(r2)^-, where r is 

the distance from the origin of the central potential. Con- 
sidering a density Ni for electrons in each quantum state 
jv^i), the magnetization of the medium can be written as 



N,q^ 



6m 



(30) 



So the coefficient am, defined by the relation M = 
"mB/^o = XmH, SO that a„i = Xm/(1 + Xm), is 



an 



6m 



(31) 



The variation of the electrons kinetic energy density 
can be written as 

«mat = l^Nd — 



E^^ 



2ms 
12m 



(32) 



Using Eq. (31), we can see that this energy density can 



be written as the second term of Eq. ([8j), 

B^ 



^mat 



2^0 



(33) 



since in this model am does not depend on the frequency 

UJ. 



VI. CLASSICAL MODEL FOR THE MATERIAL 
ENERGY FLUX AND "HIDDEN MOMENTUM" 

Here we reproduce the model of Penfield and Haus 
[25] that shows that a magnetic dipole under the action 
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of the circuit, the total linear momentum of the circuit is 



FIG. 1: Magnetic dipole in an electric field has associated 
energy flux and linear momentum. 



of an electric field can have associated energy flux and 
momentum, even if the dipole is not moving. Consider 
a rectangle of sides a and h that conducts an electric 
current / in a region of space with an electric field E. 
The situation is depicted in Fig. [T] The magnetic dipole 
moment of the object is m = ably. The electric field 
does a positive work on the positive charges that go up 
in the right side of the figure and a negative work on 
the positive charges that go down in the left side (if we 
consider that the moving charges are negative, the final 
conclusions will be the same). So a charge q moving 
to the left has an extra energy qEb with respect to the 
situation when it moves to the right. This fact generates 
an energy flux from right to left. If we have a density N of 
these dipoles, the energy flux can be found multiplying 
the energy per unity time lEh that goes to the left in 
each dipole by the length a along which the energy flows 
and by the dipole density A^. So we have 



= -lEbaNz = -E X M 



(34) 



where M is the magnetization of the medium. This result 
agrees with Eq. 

The simple model presented in this section is very dis- 
tinct from an atomic system. However, since the mag- 
netic response of an atom due to a magnetic field oscil- 
lating at optical frequencies is due to the orbital motion 
of the electrons, a similar phenomenon should occur in 
that case. 

Now let us analyze the relativistic linear momentum 
of the object. The momentum of the charges that move 
to the left in Fig. [T] is higher because they have more 
energy. The ratio between the momentum and the energy 
of a relativistic particle is v/c^, where v is the particle 
velocity. Let us consider that the current / is uniform 
along the circuit. So if we have a linear charge density 
A; with velocity vi going to the left and a linear charge 
density Xr with velocity Vr going to the right, we have 
/ = XiVi = XrVr- Calling U the energy of a particle that 
goes to the right in Fig. [ijand q the total moving charge 



XrVruU XiVia[U + qEb] 



qc^ 
m X E 



qc^ 



labE 



(35) 



So, if we have a magnetization M in a medium with 
an electric field E, there will be a "hidden momentum" 
density phid = — E x M/c^. It is worth mentioning that 
the "hidden momentum" is a relativistic effect. If we had 
considered that the momentum of a particle is mv, there 
would be no "hidden momentum" . 

In recent works, Barnett et at [5D1 [221 [30] associated 
the Abraham momentum with the electromagnetic ki- 
netic momentum and the Minkowski momentum with 
the electromagnetic canonical momentum. Considering 
the "hidden momentum" , we see that the relativistic mo- 
mentum of an atom with magnetic moment m can be 
written as Prei = Pkin — Ex m/c^, where Pkin = ™v, 
m being the mass and v the velocity of the atom. If we 
have a bunch of these atoms forming a medium, the dif- 
ference of the relativistic and the kinetic momenta of all 
atoms is equal to the difference of the total momentum 
of the field when we consider the Abraham momentum 
'Abr, and the total momentum Po.m. when we consider 



the momentum density of Eq. ([2]). So we can write 



jmcd It 

rol — ^Ahi- 



>incd I) 

kin — -t^Min 



>mcd 



(36) 



The last equality is derived in Ref. ^20j. The total ki- 
netic and canonical momenta are always conserved, so we 
conclude that, if we associate Pc.m. with the relativistic 
momentum of the field, the total relativistic momentum 
will be always conserved. 



VII. CONCLUSIONS 

In this paper we defended a particular division of the 
energy and of the momentum of electromagnetic waves in 
linear media in electromagnetic and material parts as be- 
ing the most natural. The electromagnetic parts of the 
energy density, energy fiux and momentum density are 
given by Eqs. ([T]) and ([2|, having the same form than in 
vacuum when written in terms of the macroscopic electric 
and magnetic fields. The material momentum is calcu- 
lated directly from the Lorentz force law, the material 
energy corresponds to the kinetic and potential energies 
of the charges of the medium and the material energy 
fiux is the result of the work done by the electric field on 
the charges that compose an object with magnetic dipole 
moment. 

In a previous work |23j . we showed momentum con- 
servation in various circumstances considering the elec- 
tromagnetic momentum density of Eq. ^ and calculat- 
ing the material momentum with the Lorentz force law. 
Here we showed, using reasonable models for linear me- 
dia, that we have energy conservation considering the 
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energy density of Eq. ([T]) and calculating the material 
energy density as the sum of the kinetic and potential 
energy densities of the charges of the medium. To jus- 
tify the material energy flux, we used a model introduced 
by Penfield and Haus [5^. We also argued that the mo- 
mentum density of Eq. ([2| can be associated with the 
relativistic electromagnetic momentum. It is worth to 
mention that, although we used models for linear non- 
absorptive media, we believe that the proposed division 
can be extended to any kind of media. For example, on 
an absorbing medium the absorbed energy is usually con- 
verted into heat or chemical reactions, that can also be 
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